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First observation of (anti)deuterons in deep inelastic ep scattering (DIS) measured with 
the ZEUS detector at HERA is reported. The production rate of deuterons is higher 
than that of antideuterons. However, no asymmetry in the production rate of protons 
and antiprotons was found. The (anti)deuteron yield is approximately three orders 
of magnitude smaller than that of (anti)protons, which is consistent with the world 
measurements. 



1 Introduction 

Deuterons (d) are loosely bound states whose production in high energy collisions can hardly 
be accommodated in the standard fragmentation models. The measurements of such parti- 
cles may provide an important information on the structure of fragmentation region [1] and 
on the formation of multiquark states [2]. 

In collisions involving elementary particles, several measurements have been performed in 
order to understand the production of d(d). The cross section of d in e + e _ — * qq collisions [3] 
is lower than that measured in hadronic Y(IS) and T(2S) decays, and disagrees with the 
predictions based on the LUND string model [4]. The (d) rate in e + e~ — ► qq is also lower 
than that in hadronic (pA [5], pp [6]) and in jp collisions at HERA [7], but higher than that 
in nucleus- nucleus (AA) collisions [8]. 

According to the coalescence model [1] developed for heavy-ion collisions, the d rate is 
determined by the overlap between the wave function of a proton, p, and a neutron, n, 
with the wave function of a d. The d cross section is the product of single-particle cross 
sections for p and n, with a coefficient of proportionality, Bi, reflecting the spatial size of 
the fragmentation region emitting the particles. 

Unlike studies in AA, Ap and pp, all previous measurements in collisions involving ele- 
mentary particles have been performed for d, since the reconstruction of d requires a careful 
separation of such states from particles produced in interactions of colliding beams with 
residual gas in the beam pipe and secondary interactions of particles on detector material. 
The deep inelastic scattering (DIS) events, which are characterised by the presence of a 
scattered electron in the final state, provide an ideal environment for studies of d, since the 
background from non-ep interactions is minimal. In addition, during the 1996-2000 data 
taking, the ZEUS detector had a small amount of inactive material between the interac- 
tion region and the central tracking detector (CTD). This material consisted of the central 
beam-pipe and inner wall of the CTD, with the overall thickness of 2.9 mm of Al. This leads 
to a small absorption rate of d, as well as to a small contribution from secondary deuterons. 
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Figure 1. Left: The mass spectra reconstructed using dE/dx for positive and negative 
particles. Right: The DCA distributions for: (a)-(b) particles and (c)-(d) antiparticles. 



The data sample used in the analysis corresponds to an integrated luminosity of 120 pb 1 
taken between 1996 and 2000 with the ZEUS detector at HERA. After a DIS selection, 
the average Q 2 for the data sample was about 10 GeV 2 . The particle identification was 
performed using the dE/dx measurement. The masses, shown in Fig. l(left), were calculated 
from the measured track momentum and energy loss using the Bethe-Bloch formula. 

To identify particles produced in ep collisions, the distance of closest approach (DCA) 
of the track to the beam spot in the transverse plane was used, since particles originating 
from the primary ep collisions feature small values of |Z?CM|. The DCA distribution after 
the mass cuts and an additional cut on the distance, A.Z, of the Z-component of the track 
helix to the primary vertex, is shown in Fig. 1 (right). The number of particles produced in 
ep was assessed from the DCA distribution by using a side-band background subtraction. 
The numbers of d and d after the side-band background subtraction were 177 ± 17 and 
53 ± 7, respectively. This difference was found to be unlikely related to the CTD efficiency, 
which usually leads to a larger number of negative tracks compared to positive ones: For 
example, the number of reconstructed p(p) in the data after the DCA side-band background 
subtraction was 1.52 x 10 5 (1.62 x 10 5 ). Such p— p symmetry is fully accounted for by known 
difference in the tracking efficiency for positive and negative tracks. 

Several sources of background processes for the d sample were considered [9] : events due 
to interactions of the incoming proton (or electron) beam with residual gas in the beam 
pipe (termed beam-gas interactions) and secondary interactions of particles on material 
between the interaction point and the central tracking detector. Extensive checks have been 
performed to exclude the first source. In particular, a special event selection was used for 
non-colliding electron and proton bunches. It was found that, after the DIS event selection, 
a contribution from the beam-gas events is unlikely. 

Even for a clean DIS sample, d can still be produced by secondary interactions of particles 
on inactive material. One possible source for d is the reaction N + N — > d + 7r + , where one 
nucleon, N, is produced by ep collisions, while the second one comes from inactive material 
in front of the CTD. Secondary d can also be produced in the pickup reaction p + n — » d 
by primary nucleon interacting in the surrounding material. Checks for such sources of 
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background were either negative or not conclusive due to insufficient information on the 
production cross section of the pickup process. In particular, a check was done using the 
HERAII data collected with the ZEUS detector equipped with a vertex detector. For the 
dE/dx measurement, this additional detector increased the overall material between the 
interaction point and the CTD by a factor three. As expected, the production rate of d has 
significantly increased due to a stronger contribution from spallation processes. However, 
the rate of d after the DC A background subtraction was the same as for HERAI. 
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Figure 2. Left: d/p, d/p, p/p and d/d production ratios as a function of pr/M. Right: 
Comparison of the B2 values extracted from DIS with other world measurements. 



The detector-corrected production ratios as a function of pt/M are shown in Fig. 2 (left). 
For the antiparticle ratio, there is a good agreement with the HI published data for photo- 
production [7] , as well as with pp data [6] . The production rate of d is higher than that of d 
at low pt- If (anti)deuterons are produced as a result of a coalescence of two (anti)nucleons, 
then one should expect that the d/d ratio is approximately equal to the (p/p) 2 ratio for the 
same pt per (anti)nucleon, assuming the same source radius for particles and antiparticles. 
In this case, one should expect d/d ~ 1, since the measured p/p ratio is consistent with 
unity. Under the assumption that secondary interactions do not produce an enhancement 
at DC A = for the d case, the result would indicate that the relation between d/d and 
(p/p) 2 expected from the coalescence model with the same B2 for particles and antiparti- 
cles does not hold in the central fragmentation region of ep DIS collisions. In terms of the 
coalescence model, the d production volume in momentum space is larger than that for d. 

For collisions involving incoming baryon beams, there are several models [10] which 
predict p— p production asymmetry in the central rapidity region. The predicted asymmetry 
can be as high as 7% due to the presence of the incoming proton. As shown in Fig. 2 (left), 
the experimental data for p(p) are not sufficiently precise to confirm such expectations. 

There are no predictions for the d — d asymmetry. It is possible that theoretical expec- 
tations for such compound states are different than those for the p — p asymmetry, since 
d(d) are not contaminated by a large contribution from the standard baryons produced in 
quarks and gluon fragmentation. 

The production of d(d) was studied in terms of the coalescence model. The B2 values 
extracted in the region 0.3 < pr/M < 0.7 are shown in Figure 2(right). While the mea- 
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surement of B2 has larger experimental errors than those in the studies of the production 
asymmetries, it is still seen that B2 for d tends to be higher than that for d. The values 
of B2 for d are in agreement with the measurements in photoproduction [7], but disagree 
with the B2 measured in e + e~ annihilation [3] at the Z resonance. In contrast to heavy-ion 
collisions, where B2 strongly increases as a function of pt due to an expanding source, the 
B2 measured in ep does not strongly depend on pr/M [9], 

In summary, the first observation of (anti)dcuterons in ep collisions in the DIS regime is 
presented. The yield of d(d) is three orders of magnitude smaller than that of p(p), which 
is in broad agreement with other experiments. The production rate of p is consistent with 
that of p, while the production rate of d is higher than that for d for the same kinematic 
region. 
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